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In this paper, we study the vanishing viscosity hmit of one-dimensional isentropic 
compressible Navier-Stokes equations with density-dependent viscosity, to the isen- 
tropic compressible Euler equations. Based on several new uniform estimates to 
the viscous systems, in addition to the framework recently established by G. Chen 
and M. Perepelitsa |10] . we justify that the finite energy solution of the isentropic 
compressible Euler equations for a large class of initial data can be obtained as the 
inviscid limit of the compressible Navier-Stokes equations even when the viscosity 
depends on the density. 
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1 Introduction 

When the fluid density experiences large scale dropping, especially when vacuum is con- 
cerned, the motion of isentropic compressible viscous fluids is modeled by the following 
compressible Navier-Stokes equations with the density-dependent viscosity, in the Eule- 
rian coordinates. 
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where and denote the density and the velocity of the fluid, respectively, nf = p^u'^ 
represents the momentum, p = p{p) is pressure function of the density. In this paper, we 
consider the polytropic perfect gas, i.e. 

p{p) = np^, 

where 7 > 1 is the adiabatic exponent, and the constant k is chosen as k = ~^ up 
to a scaling. While e > is adpated to the system as the controlling parameter on the 
amplitude of viscosity, for which we assume e G (0,£o] for some fixed Eq > without loss 
of generality; a > is a constant which models the dependence of viscosity on density. 

When a = 0, (1.1) reduces to the classical compressible Navier-Stokes equation, called 
CNS. The case of a > occurs for non-uniform gases [5], and (1.1) can be formally 
derived by Chapman-Enskog expansion from the Boltzmann equation for (at least) hard 
sphere model and cut-off inverse power force model. A formal derivation can be found in 
[39] . It is also interesting to note that when a = 1 and 7 = 2, (1.1) recovers the "viscous 
Saint- Venant" system for shallow water without bottom friction |T9], see also [11]. In this 
paper, we will focus on (1.1) with positive a, for which we call it a-CNS, distinguishing 
from the case of a = 0, which is called CNS. On the other hand, the studies in [23] and 
[37] indicate the failure of CNS at vacuum and the validity of a-CNS at least at the level 
of local well-posedness theory. We therefore devote our efforts to this model in current 
paper. 

We now consider the Cauchy problem of (1.1) when the far fields of the fluid are away 
from vacuum. Namely, we shall study the a-CNS (1.1) with the following initial data 

p'{0,x) = pl{x) > 0, u%0,x) = u^ix), (1.2) 

such that 

lim {pq{x) , Uq{x)) = {p^,u^), with p^ > 0. 
x— ^±00 

In the past decades, the study of the mathematical theory on (1.1)-(1.2) has attracted 
a lot attention. Many interesting results were established for the local and global existence 
of both classical and weak solutions, we refer the readers to some of them such as, [2], [3], 
[E], [18], [22], [27], [28], [29], [30], [33], [37], [39], [30], [31] and [51]. It is equally interesting 
to study the inviscid limit for (1.1)-(1.2) as e — )■ toward the following one-dimensional 
isentropic Euler equations 

Pt + {pu)x = 0, 

{pu)t + {PU"^ +P{p))r, = 

It is a general belief that the physcial weak solution of (1.3) can be obtained in such a 
process, see [3], where a vanishing artificial viscosity limit for general hyperbolic system 
with small BV data is proved. This problem is closely related to the existence of weak 
solutions to (1.3) through a limitting process of physical approximation. In this paper, 
we will address this problem and study the vanishing viscosity limit for (1.1)-(1.2). 
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(1.3) 



In BV framework, when the initial data is away from vacuum, the existence of global 
BV solution to (1.3) was established by [H] for 7 > 1 and by [15] for 7 = 1 using Glimm's 
method. (1.3) shows singular behavior when vacuum occurs which causes difficulties to 
mathematical analysis. It is still a major open problem on how to perform BV estimate 
when the solution may contain vacuum states. Instead, the L°° framework is successfully 
achieved using the theory of compensated compactness [13], [S2]. The existence of L°° 
weak entropy solution of (1.3) was established by [I2] for 7 = 1 + n > 2; by [13] for 
7 G (1, §]; by [32] and [SB] for 7 > |; and finally by [2B] for 7 = 1. Recently, [32] further 
constructed the finite-energy solutions to the isentropic Euler equations with finite-energy 
initial data. We remark that these results are achieved through the vanishing artificial 
viscosity. 

The problem of vanishing physical viscosity limit is more subtle and the progress 
has been less satisfactory, and the problem of vanishing viscosity limit of Navier-Stokes 
equations to Euler equations has been open for long time, though some interesting results 
are proved when restrictive initial data is assigned, see [23] and [53]. Recently, G. Chen and 
M. Perepelitsa [10] proved that the solutions of Navier-Stokes {a = 0), whose viscosity 
is independent of density, converge to the finite energy solution of Euler equations as 
viscosity vanishes. This is a major breakthrough in this aspect. 

Inspired by [10], we study the problem of vanishing viscosity for the a-CNS (1.1)-(1.2) 
in this paper with positive a. It is clear that for any fixed positive e, the visocity coefficient 
with positive a experiences degeneracy near vacuum states. An obvious obstacle is the 
dissipation term in the energy identity contains only the weighted norm of velocity gradi- 
ent which degenerates at vacuum. Such a singular behavior causes the major difficulty in 
the analysis and introduced the different behavior of solutions compared with CNS where 
a = 0. The analysis exibits quite different fiavor and requires very different ingredients. 
Fortunately, by a deep observation, we obtained several key uniform estimates. Based 
on these uniform estimates and the framework of [TU], we are able to show that, when 
viscosity parameter £ tends to zero, the solutions of a-CNS fll.ip -(1.2) converge to the 
finite-energy solution of Euler equations for general initial data. 

We now prepare to state our main result. 

A pair of functions {ri{p,u),q{p,u)), or {r]{p,m),q{p,m)) for m = pu, is called an 
entropy-entropy fiux pair of system (1.3), if the following holds 

[r]{p,u)]t+ [q{p,u)]^ = 0, 

for any smooth solutions of (1.3). Furthermore, rj{p,m) is called a weak entropy if 

r](0, u) = 0, for any fixed u. 

An entropy r]{p, m) is convex if the Hessian V'^r]{p, m) is nonnegative definite in the 
region under consideration. 
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From [36], it is well known that any week entropy {rj, q) can be represented by 



(1.4) 

q^{p, pu) = q^{p, m) = / (6's + (1 - 9) u) x{p;s-u)ij (s) ds. 

Jr 

where the kernel is x(p; s — u) = [p^^ ^ {s — ufY+i ^ ~ "^^-i) ~l' ^^'^ ^ ~ ^T^- ■^^'^ 
instance, when -0(5) = |s^, the entropy pair is the mechanical energy and the associated 
flux 

2 

Tfi TYl 

r^*[p^rn) = — + e{p), q*{p,m) = — + me'{p), (1.5) 

where e(p) = t~[P"' represents the gas internal energy in physics. 

Let {p{x) ,u{x)) be a pair of smooth monotone functions satisfying (p(x) ,u{x)) = 
{p^,u^), when ±x > Lq for some large Lq > 0. The total mechanical energy for (11. ip in 
M with respect to the pair of reference function [p (x) , u (x)) is 

E[p,u]{t)= / {r]* {p,m) — rj* {p,m) — Vrj*{p,fh) ■ {p — p,m — m)) dx, (1.6) 
Jr 

where rh = pu. After some calculations, we obtain that 

E[p,u](t)= I ( -p(t,x)\u(t,x) — u{x)f + e* {p{t,x),p{x))] dx (1.7) 



where e* (p, p) = e(p) — e(p) — e' (p) {p — p) > 0. 

Definition 1.1 Lei (po,Uo) &e g'zwen initial data with finite- energy with respect to the end 
states {p'^jU'^) at infinity, and E[pQ,uo] < Eq < oo. A pair of measurable functions 
{p,u) : — is called a finite-energy entropy solution of the Cauchy problem (1.3) if 
the following holds: 

(i) The total energy in bounded in time: There is a bounded function C {E,t), defined 
on ]R+ X and continuous in t for each E G M^, such that, for a.e. t > 0, 

E[p,u] {t)<C{Eo,t)- 

(a) The entropy inequality: 

ip,u\ + q^ {p,u)^<0, 

is satisfied in the sense of distributions for all test functions ip (s) G {±1, ±s, s^}; 
(Hi) The initial data (po,tto) are attained in the sense of distributions. 

We now state our main conditions on the initial data (1.2), which is motivated from 
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Condition 1.1 Let {p{x) ,u{x)) be some pair of smooth monotone functions satisfying 
{p{x),u{x)) = (p^ (x) , (x)) when ±x > Lq for some large Lq > 0. For positive 
constants Eq, Ei and Mq independent of e, and Cq > the initial functions {Po,Uq) are 
smooth satisfying the following properties 

(i) Po > Cq > 0, / Po (x) \uq (x) — u{x)\dx < Mq < oo ; 
Jr 

(a) The total mechanical energy with respect to {p,u) is finite: 

j Qpo ko - i.^)? + (Po i.^) > P i.^))^ dx <Eq < oo; 

(^'^^) / ^-rfr^dx <Ei< oo; 
Jr Po [x) 

(iv){pl (x) , Pq (x) (x)) — > (po (x) , Po (x) Uq (x)) in the sense of distributions as e ^ 
0, with Pq (x) > a.e.. 

Our main results are stated in tlie following Theorem. 

Theorem 1.1 Assume |<a<7, 7>1. Let {p'^,u'^), rrf = p'^u'^ be the solution of 
the Cauchy problem M.l\] -(1.2) with initial data {pq,Uq) which satisfies Condition 1.1 for 
each fixed e > 0. Then, when £ — )■ 0, there exists a subsequence of (p"^, rrf) that converges 
almost everywhere to a finite-energy entropy solution (p, m) to the Cauchy problem 
with initial data {po,Po'^o) for the isentropic Euler equations. 

Remark 1.1 Due to some technical difficulty, we can only prove the result for | < a < 7. 



In fact, for many physical gases the Chapman- Enskog viscosity predicts that a > see 
J^, and 13^ . Our condition | < a < 7 «s valid for many physical cases including the 
shallow water model, but it did not cover the case of monoatomic gas where 7 = | and 
a = |. It is very interesting to prove the result for a G [0,|], which will be addressed 
later. 



One important basis of our proof for Theorem 11.11 is the following compactness theorem, 
established in 1101. 



Theorem 1.2 (Chen-Perepelitsa [ID]) Let G C^{R), (77^, q"^) be a weak entropy 
pair generated by i/j. Assume that the sequences (p^(x, t), m^(x, t)) defined on M x with 
rrf = p'^u^ , satisfies the following conditions: 

(i). For any —00 < a < 6 < 00 and all t > 0, it holds that 

[ [ {p'y+^dxdr <Cit,a,b), (1.8) 

Jo Ja 

where C{t) > is independent of e. 

(a). For any compact set K C M, it holds that 

[ [ {py^^ + p'\u'fdxdT <C{t,K), (1.9) 
Jo Jk 
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where C = C{t,K) > is independent of e. 

(Hi). The sequence of entropy dissipation measures 

rf^p"^, rrf)t + g'''(p^, rrf)^ are compact in _ff^~^(M^). (1-10) 

Then there is a subsequence of {p'^ ,rrf)( still denoted as {p'^^rrf)) and a pair of measurable 
functions (p, m) such that 

(p'^ ,m'^) {p,m), a.e. as e 0. (l-H) 

In section 2 and section 3 below, we will verify conditions (11. Sp . (II. 9p and (ll.lOp to 
prove our main theorem 11.11 

The rest of this paper is arranged as follows. In section 2, we make some new uniform 
estimates for the solutions of Navier-Stokes equations (11.11) which are independent of e. 
These estimates are essential to show the convergence of the vanishing viscosity limit to 
the Euler equations. In section 3, using the estimate we obtained in section 2, we prove the 
i/;~^(R^)— compactness for the solutions of ( 11. ip . In section 4, based on the framework 
in [To], we prove our main Theoren ll.il 

2 Uniform Estimates for the Solutions of a-CNS 

First, we assume that (p^,w^) is the global smooth solutions of Navier-Stokes equations 
dLlD-dLSD, satifying 

p%x,t) > c%t), for some c%t) > (2.1) 

and 

hm (p^^x^)(x,t) = (p±,M±). (2.2) 

X— >±oo 

For the existence of global smooth solutions, the reders are referred to [28], [39] and |41j . 
Based on the above preparation, we now make some new the uniform estimates with 
respect to e for the solutions (p^,M^) of the a-CNS (ll.ip -( fL2l) . 

For simplicity, throughout this section, we denote (p, tt) = (p^,M^) without causing 
confusion and C > denote the constant independent of e. 

Lemma 2.1 (Energy Estimates) Suppose that < a < 7, and E[po,uo] < Eq < 00 for 
some Eq > independent of e. It holds that 

sup E[p,u]{T) + e [ [ p'^uldxdr <C{t), (2.3) 

0<r<i Jo Jr 

where C{t) depends on Eq, t, p, and u, but not on e. 
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Proof. From the definition, we have 

^^^^^ - ^ /" r]*{p,m) - r]*{p,m) dx - I Vr]*{p,m){pt,mt) dx. 



dt dt . 
Since {'r]*,q*) is an entropy pair, we have 

r]*{p, m)t + q*{p, m)^ - er]*^{p, m){p°'u^)^ = 0. 
Integrate (12 .Sp with respect to x over M, we obtain 

r]*{p, m) — ri*{p, fh) dx + e / p^u^ dx = q*{p~ ,m~) — q*{p^ , m 
: Jr 

Since we have 

e*(p,p)>p(p^-pT, 9 



d_ 

di 



.e^2 n 7-1 



2 

Utilizing (12. 7p . we obtain 

Vr]*{p,m){pt,mt)dx\ 

Vr]*{p, rh){m^, {p{p) + pu^ - ep"'u^)^)dx\ 
Vr]*{p, fn).j,{m,p{p) + pu^ - ep^Ur,) dx\ 



< I / p'^ul dx + C I p\u- dx + C I {p + p{p) + p'')dx + C 



Lo 



<C + CE + ^ I p'^uldx, 



where we have used 

"Lo 



(p + p{p) + p'^)dx < CE, for < a < 7. 

-Lo 

Substituting and ([23D into ([MD, we obtain 

■ ^ / P^'ul dx<C + CE, 



dt 4 

Then Gronwall's inequality implies Lemma l2m 

Remark 2.1 Since vacuum could occur in our solution, the inequality 



t 

' p'^ul dxdr < C{t) 
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in (12 .Sp is much weaker than the corresponding one 

u^dxclr < C{t). 







in fT^ - This will cause a great difficulty to prove Lemma WT^ below, which is an essential 
step to verify the condition i) of Theorem M.^ i.e. (11.81) . 

We now derive some higher order estimates. 
Lemma 2.2 IfO<a<^, and {po{x) , uo{x)) satisfies 

2 



\pox{x)\ 



^2 I iruxv-.i d^^E,< CX), 



Po[x) 

for some Ei independent of e. Then, for any t > 0, it holds that 

^2 I p2a-3p2^ dx + e [ [ p'^^^-^pI dxdr < C{t), (2.10) 
Jr Jo Jr 

where C{t) depends on Eq, Ei, t, p, u, hut not on e. 
Proof. Through (11.11) . we have 

Pxt + PxxU + 2p^Ux + pu^x = 0, 
put + puu^ + p{p)^ = ep°u^^ + eap°'~^u^pr,. 

Multiplying (I2.1ip i with p^"~^p^., after some calculation, we obtain 



^2 2a-3 2\ /^2 2a-3 2\ 

e_p p^^ ^(._p up^^ ^^^^«+7-3^2 



+ ( ep"-VxW ) - ( sp'^-'ptu ] = ep'^ul 



Integrating (I2.14p over [0,t] x M, we obtain 

^2 f p!!^ + eK f I p'^^^-^pI dxdr 







< C{t) -e [ [ (p'^-^p^u) dxdr. 

Jo Jr\ J t 



(2.11) 



{^^)+ (^^^) + ap'^-'plu. + p'-'p.u.. = 0. (2.12) 

From (I2.IIP 9 X p^^^p^,, after some calculation, we reach 

p"~VxMt + p'^^^up^u^ + kp'^^^-'^pI = eap'^"-^plu^ + ep^^-Vx^xx- (2.13) 
The combination e^l^J^ + g (l2J[3D gives that 



(2.14) 



^2 / Po pIx ^^^^ I f p^uldxdr-e [ I (p"~Vx^) dxdr (2.15) 



Noticing that 



p"-Vx« ) = ( P""Vx(^ -u)] + - ( ) , (2.16) 



t \ / t ^ \ / t 



e 



{p''),u dx = — p'^ujx + - (p+)"w+ - {p^Tu 

Li 



JO 



(2.17) 



Integrating f l2.16p with respect to x over M, we have 

el i p^'^pxU J dxcir 
Jo Jr\ ) t 

= e p°'~'^pnc{u -u) dx - e I Po"Vox(mo - u) dx 
Jr Jr 

+ - {p'^)xudx + - {Po)xU dx 
<Cit), 

where we have used the following estimates fl2.18l) - fl2.20l) 

e / p°''^^ px{u — u) dx 
Jr 

< — [ p'^'^pldx + C I p{u~ uf dx (2.18) 
8 Jr Jr 

° Jr 

e I pVpUuo -u)dx<^ [ pI'^-'pIJx + C(t), (2.19) 
Jr ° Jr 



< — I p^dx + C{t) 

Lo 

< C(t), (2.20) 



Substituting (l2T71) -( l2:20ll into ( ]2J[5l) . we obtain LemmeEll 

The following higher order integrability estimate is crucial in compactness argument. 

Lemma 2.3 // the conditions of Lemma \2. 1\ hold and < a < 7, then for any —00 < 
a < 6 < 00 and all t > 0, it holds that 

[ [ p''+^ dxdr <C{t,a,b), (2.21) 

Jo J a 

where C{t) > depends on Eq, a, b, 7, t, p, u, hut not on e. 
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Proof. Choose 

w{x) G C^(M), < w{x) < 1, w{x) = 1 for X G [a, b], and supp{w} = (a — 1, + 1) 
By fll.iP n X w, we have 

{p{p)w)x = -{pu^w)x + + - {pu)tw + e{p°'Uxw)x - ep^u^w^ (2.22) 

Integrating (12.221) with respect to spatial variable over (— oo,x), we obtain 

p{p)w = —pu^w + ep^u^w ~ (^J P'^'^ 

+ [{pu^ +p{p))wx- ep^u^Wx] dy. 



Multiplying (I2.23P by pw, we have 
pp{p)w^ = —p^u^w"^ + ep°'~^^Uxw'^ — { pw I puw dy 

\ J — oo / 

- {pu)^w I puw dy + pw I [{pu^ + p{p))wx - ep'^UxWx] dx 
= ep°'~^^UxW^ — { pw I puw dy \ — ( puw I puw dy 

\ J -co J t V ^— oo J 

/X rx 
puw dy + pw / [{pu^ + p{p))wx - ep^'uxWx] dy 
-oo J — oo 

Integrating (12.241) over (0,t) x M, we have 
np^^^w"^ dxdr 

Pqw I pouowdy 







oo 

X 



jdx + i puwx / puw dy \ dxdr 

J Jo Jr\ J~oo J 



pw I [{pu'^ +p{p))wx - ep°'UxWx]dy ) dxdr 

oo 



(2.23) 



(2.24) 



(2.25) 



Let 

A = {x : p{x,t) > p}, where p = 2 max{p+, p— }, (2.26) 
then we have the following estimates by (12.31) 
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By (I2.26p . we know that for any there exists a point xq = Xo{x,t) such that 

— a;o| < d{t) and p{xo, t) = p. Here we choose (5 = a + > 0, 

sup ep^ {x,t) < ep^ + sup ep^{x,t) 

xGsupp{w} x^supp{w}nA 

<2ep^ + sup \ep^{x,t) - ep'^{xQ,t)\ 

x(isupp{w}f\A 

fXo+d{t) 

I .,nr. / IRW^n^^h 



<2£p^+ sup / IPWep^^'p^l dx 

xGsupp{w}nA J xo—d{t) 
fb+l+2d(t) 

<2ep^+ mep^-'p,\dx (2.28) 



'a-l-2d{t) 
''b+l+2d{t) 

< 2epf + / l/^lp^/'-^^+i dx + I e^p^'^-'pl dx 

Ja-l~2d{t) Jr 
l'b+l+2d(t) 

< C{t) + / p^ dx 

Ja-l-2d{t) 

<C{t) 

Using (12.28p . the first term on the right hand side of (12.25P can be estimated as 
following 

/ p°'^^Uxn?dxdT 
Jr 

<£ [ I p^'+^w^ dxdT + e I I p'^ul dxdr 
Jo Jr Jo Jr 

< C{t) +e [ I p'^^'^n? dxdr 
Jo Jr 

/ / p'^w'^ dxdr, if 2a + 2 < ^ 
J,o Jr 

sup £p^(x,r) f p^^^^-f^w^ dxdr, if 2a + 2 > ^9 
?uppiw} Jr 



< 




(2.29) 



x&supp{w} 

t 



< C{t) + C{t) I I p2"+2-/3^2 ^^^^ 

Jo Jr 

< C{t) + 5 [ I p^+^w^ dxdr. 

Jo Jr 

Here we have used the fact 2a + 2 — /? < 7 + 1 for 7 > 1. 
By Lemma 12.11 and the Holder inequality, we obtain 



puwdy\ < / \pu\dy 

J supp{w} 



1 1 

2 / r \ 2 



pdy\ I pu'dy] <C{t). (2.30) 

SUpp{w} J \J SMppjio} 
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Then it follows that 

pw I puw dy ) dx\ + | / ( PoW I PqUqW dy ) dx\ (2-31) 



/ / I puWx / puw dy j dxdr] < C{t). (2.32) 

Jo Jr \ J-oo / 



Similarly, we have 

ft /• / fX 



and 



pw / {pu^+p{p))w^ dy ] dxdrl < C{t), (2.33) 

V J-oo 



\e i pw / p'^UxWx dy j dxdT\ 

Jo Jr V J-oo J 

< e / / I pw / p^l'Uj;! dy ] dxdr 







(2.34) 

<e I ( f pwdxV f p-ul dy+ j p^wl dy\ dr 
<C{t) 

Substituting fl2:29|) . fl23TD - fl23D into (12:251) and noticing the smallness of 6, we 
proved Lemm£(231 

Lemma 2.4 Suppost that | < « < 7, and (po(3^)) ''^0(2^)) satisfy the conditions in the 
Lemmas \2.1\ \2.Bt \2.3\ . Furthermore, assume that for some Mq > independent of e, the 
following 



Po{x) \uq{x) — u{x) \ dx < Mq < 00, (2.35) 
holds. Then for any compact set K C M, it holds that 

[ [ p^+^ + p\uf dxdr <C{t,K), (2.36) 
Jo Jk 

where C = C{Eo, Ei, Mq, t, K) > independent of e. 

Proof. First, we introduce a useful result about the entropy pair, see [36] for details. 
Taking ip^iw) = ^w\w\, then there exists a positive constant C > 0, depending only on 
7 > 1, such that the corresponding entropy pair {ri*,q*) = {rj^* ,q^*) satisfies 

(\v%P,u)\ < {p\u\^ + p^), 

^ ' - ' ' for all p > and u E M. 2.37 

\v:^ip,u)\<C{\u\+p'), 

J^mm(P)«)l < Cp~\ 



12 



If r)^ is regarded as the function of (p, u), we have 

\vUp,^)\<C, forallp>Oand«GM. 

\vt,,ip,u)\<Cp'~\ 

For this weak entropy pair [r]*, q*), we note that 

v'iP, 0) = v'M 0) = 0' i'(p^ 0) = ^oP''^' I l^ni - ^']+ 

and 

r7;;(p,0)=/3p^ with/3:= /" ds. 
Taylor expansion imphes 

rf'i^p, m) = Pp^m + r(p, m), 

with 

r(p, m) < Cpu^ . 

for some constant C > 0. Now we introduce a new entropy pair (r), g), 

7}(p, m) = rj*{p^ m — pu~), q{p, m) = q*{p, m — pu~) + u~ri*{p, m — pu~ 
which satisfies 

fi{p, m) = j3p^^^{u — u~) + r(p, p{u — u~)), 
r{p, p{u — u^)) < Cp{u — u~Y. 

Integrating (11. IP -^ + (I1-1P 2 x Vm over (0, t) x (— oo, x), we obtain 

(?7(p, m) - ?7(po, mo)) dy + I q*{p, p{u - u)) + u~r]*{p, p{u - u)) dr 

3 Jo 

= tq*{p~,0) + e fj^p'^u^ dr - e / {fjmuP'^ul + fimpP"pxUx) dydr. 

Jo Jo J -oo 

Utilizing f l2.38p . we have the following estimates 

fimuP''uldydT\ <Ce f [ p'^uldydr < C{t), 
Jo Jr 

pt px ft f 

I I Vmpp"^ pxUxdydT\ <Ce / p^^VlPx^^l dydr 
Jo J-oo Jo Jr 

<Ce [ [ p^uldydT + Ce [ [ p^'^^'^pl dydr < C{t). 




'0 J-oo 
ft nx 
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Substituting f l2.43p and f l2.44p into (12 .42 p . then integrating the result with respect to 
X over K and using ( I2.37p . we obtain 



pa+i _j_ — u \ dxdr 



JK 

< C{t) + c 

t 

+ Ce 



t 

JK 

a+e\ 



\ri*{p, p{u — u))\ dxdr + Ce / / p°'\ux\\u\ dxdr 



p \ux\ dxdr +2 sup 
'0 JK Te[o,t] 

Applying Lemma 12.11 it is easy to get 



JK 

Vipiy,^), {pu){y,T)) dy ] dx 



(2.45) 



K \J-oo 



\r]*{p,p{u - u))\ dxdr < C{t). 



(2.46) 



JK 



Now Cauchy-Schwartz inequality and (I2.28P lead to 



£ / / P^'^^lu^l dxdr < Ce I I p'^u'^ dxdr + Ce I I p'^^'"' dxdr 



JK 



a„.2 



JK 



JK 



< C{t) + C{t) 

< C{t). 

Noticing that 3/3 > 3a - 2 and fl2:28|) . we have 



7-1 



p 2 dxdr 



JK 



(2.47) 



JK 



p°'\ux\\u\ dxdr < -e 

2 Jo 



I ft f 1 /"* 

p'^ul dxdr + -e I 
K 2 Jo 



p°'u^ dxdr 



K 



< C{t) + C{6)e^ [ [ p=^"-2 dxdr + 6 [ [ p\uf dxdr 

Jo J K Jo J K 

< C{t) + [ I p^^ dxdr + 5 [ [ p\u\^ dxdr 

Jo J K 

< C{t) + 6 [ [ p\u\^ dxdr 



JK 



(2.48) 



JK 



where S is small enough which will be determined later. 

Now we estimate the last term on the right hand side of fl2.45p . (11. ip implies that 



{pu - pu )t + {pu^ + p{p) - puu )^ = e{p°'u^). 
Integrating (I2.49P over [0,t] x (—00,2;) for x G we obtain 



(2.49) 



p{u — u ) dy 



Pq{uq~u ) dy - {pu^ + p{p) - puu -p{p )) dr 



+ e i p^Ux dr. 



(2.50) 
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< 



< 



-oo 

X 



On the other hand, 

v{p{y,T): {pu){y,r)) dy 

{flip, pu) - I3p^^^{u - u)) dy 
r{p,p{u - dy 

) 

+ P{p~f I p{u-u-)dy 

J — oo 

<C{t)+P{p-y r p{u-u-)dy 

J — oo 

which, together with fl2.35p . Lemma l2m -Lemm8 [2.3l and fl2.50p . imphes that 



pp'^+^{u-u-) dy 

) 

Pi^p^ -(^p^Y)p{u-u^) dy 



(2.51) 



K 



fi{p{y,T),m{y,T))dy 



dx < C{t). 



(2.52) 



Now, if one chooses 5 small enough, then substitutes (I232D . ([MSD, (12^48]) and fl2^ 
into (12.451) . the proof of Lemma 2.4 follows. 

Remark 2.2 In the uniform estimates above, we have required that | < a < 7, and the 
initial functions (pq(x), Mq(x)) satisfy 

(i) Pq{x) > 0, / Pq{x)\uq{x) — u{x)\dx < Mq < 00; 
Jr 

(a) The total mechanical energy with respect to {p,u) is finite: 
1 2 

-pl (x) \ul{x) - u{x)\ + e* {pI{x),p{x)) dx =: Eq < oo; 



(Hi) 



^-%^.dx <E,<oo. 
Po 



(iv) Mq, Eq, El are independent of e. These conditions are essential parts of Condi- 
tion 1 in section 1. We remark here that the limit of the functions satisfying the conditions 
i)-iv) is very general, including a wide class of L°° functions with finite energy and may 
contain vacuum. It is obvious that the above limit can serve as the initial data of isentropic 
gas dynamics for the existence of finite energy solutions. We also note that the condition 
(Hi) is slightly weaker than the corresponding one in fW^ near vacuum. We refer to fT^ 
for further details. 



3 H~l(Rl)- Compactness 



lac 

In this section we will use the uniform estimates obtained in the previous section to prove 
the following key Lemma, which states the if^~^(M^)— compactness of the approximate 
solution sequence. 
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Lemma 3.1 Let | < a < 7, "^/^ G Cq (M) , {r]^, q^) be a weak entropy pair generated by ifj. 
Then for the solutions withrrf = p^u^ of Navier- Stokes equations (1.1)-(1.2), the 

following sequence 

r]'''{p'', nf )t + nf)^ are compact m if;^](M^) (3.1) 

Proof. In order to prove this lemma, we first introduce the following results for the 
entropy pair {vj^.q^) generated hy ip E Co(M), and see [10] for details. 

For a function : M — M, compactly supported on the interval [a, 5], we have 

supp{Tf}, supp{q^} C {(/9, m) = (p, pu) : u + p^ > a, u — p^ <b^ . (3.2) 

Furthermore, there exists a constant C.^ > such that, for any p > and -u G M, we have 

(i) For 7 G (1,3], 

\v^{p,m)\ + \q^{p,m)\<C^,p. (3.3) 

(ii) For 7 G (3, +00), 

\vHp,^)\ < c^p, \qHp.^)\ < c^P + p'"-'). (3.4) 

(iii) If rj'^ is considered as a function of (p, m), m = pu, then 

l^m(P>"^)l + \pvtmiP,^)\ < (3-5) 

and, if 77^ is considered function of {p,u), then 

\vLip,pu)\ + < c^- (3.6) 

Now we are going to prove the lemma. 

A direct computation on (11.11) -^ x r]^[p'^,rrf) + fll.iP n x rlf^^p"^ ,rrf) gives 

r/^(p^ m% + g^/, "^'). = ^(^:^;(P^ - ^ryL(P% m^)(/)"«)^ 

-^r;:(;,(p^m^)(p^)v^.< (3.7) 

Let K C M be compact, using (13.61) and Cauchy-Schwartz inequality, we have 



ritu{p'MWm<f + \ritM^^')^PTpl<\ dxdt 
<Ce [ [ [p'Yiulf dxdr + Ce [ [ {p'Y+^-^{plf dxdr ^^'^^ 

Jo JK Jo Jk 

<C{t). 
This implies that 

- er^K/, rrf)ipTiuiy - evt/p^ m%pT pK are bounded in L\[Q, T]xK), (3.9) 
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and thus it is compact in W^^^'^^ {M."^) , for 1 < pi < 2. 
Moreover, noticing \ri^{p^ , p^u^)\ < C^, we have 



t ■ 



JK 



er]l{p',m'){p')''ul] dxdt 



< e^{p'')^\ul\^dxdt 

Jo J K 



< 



Cel [ [ {p'yiull^dxdt + Cel [ [ {p')^"dxdt ^^'^^^ 
Jo Jk Jo Jk 

< C{T, K)e^ + Cel [ [ (py^^dxdt 

Jo Jk 

< C{T,K)e^ -> as £ ^ 0. 
Then([3lOD and ([33]) yield that 

rf{p'',rrf)t + q'^'ip" ^m")^ are compact in Wi~^'^'^(M.\) for some 1 < p2 < 2. (3.11) 

On the other hand, using the estimates in fl3.3l) - fl3.4p and Lemma l2.H -Lemma 12.4^ we 
have 

'rf{p'^,nf),q^{p'^,'rrf) are uniformly bounded in m^{R\) for p3 > 2, (3.12) 
where ps = 7 + 1 > 2 when 7 G (1,3]; and = > 2 when 7 > 3. This yields that, 

r7^(p^ m")t + q'^ip", m")^ are uniformly bounded in W'J'^^'iRl). (3.13) 
Then (13. lip and (13.130 implies Lemma [3.11 

4 Proof of Theorem 

Proof of Theorem 11.11 From Lemmas I2.1H2.4I and the compactness estimate Lemma 
13. H we have verified the conditions (i)-(iii) of Theorem 11.21 for the sequence of solutions 
(p^, m^). Basing on Theorem II. 21 there is a subsequence (p^, m^)(still denoted as (p^, m^)) 
and a pair of measurable functions (p, m) such that 

{p^^m") {p,m), a.ee-^0. (4.1) 

It is easy to check that (p, m) is a finite-energy entropy solution (p, m) to the Cauchy 
problem (II. 3p with initial data {po,PoUo) for the isentropic Euler equations with 7 > 1. 
Therefore, the proof of Theorem 11.11 is completed. 
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